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ABSTRACT 

Context. Star forming regions may share many characteristics, but the specific interplay between gravity, magnetic fields, large-scale 
dynamics, and protostellar feedback will have an impact on the star formation history of each region. The importance of feedback 
from outflows is a particular subject to debate, as we are yet to understand the details of their impact on clouds and star formation. 
Aims. The Pipe Nebula is a nearby molecular cloud hosting the B59 region as its only active star-forming clump. This paper focuses 
on the global dynamics of B59, its temperature structure, and its outflowing gas, with the goal of revealing the local and global impact 
of the protostellar outflows. 

Methods. Using HARP at the JCMT, we have mapped the B59 region in the / = 3 ^ 2 transition of ^^CO to study the kinematics 
and energetics of the outflows, and the same transitions of ^^CO and C^^O to study the overall dynamics of the ambient cloud, the 
physical properties of the gas, and the hierarchical structure of the region. 

Results. The B59 region has a total of ~ 30 M© of cold and quiescent material, mostly gravitationally bound, with narrow line widths 
throughout. Such low levels of turbulence in the non- star-forming regions within B59 are indicative of the intrinsic initial conditions 
of the cloud. On the other hand, close to the protostars the impact of the outflows is observed as a localised increase of both C^^O line 
widths from ~ 0.3 kms~^ to ~ 1 kms~\ and ^^CO excitation temperatures by ~ 2 - 3K. The impact of the outflows is also evident 
in the low column density material which shows signs of being shaped by the outflow bow shocks as they pierce their way out of 
the cloud. Much of this structure is readily apparent in a dendrogram analysis of the cloud and demonstrates that when decomposing 
clouds using such techniques a careful interpretation of the results is needed. 

Conclusions. The low mass of B59 together with its intrinsically quiescent gas and small number of protostars, allows the identifica- 
tion of specific regions where the outflows from the embedded sources interact the dense gas. Our study suggests that outflows are an 
important mechanism for injecting and sustaining supersonic turbulence at sub-parsec size scales. We find that less than half of the 
outflow energy is deposited as turbulent energy of the gas, however this turbulent energy is sufficient to slow down the collapse of the 
region. 

Key words. Stars: formation, protostars; ISM: clouds, jets and outflows, kinematics and dynamics, individual objects: B59 



1. Introduction 

As stellar nurseries, molecular clouds contain the ingredients for 
star formation. However, the physical processes which prevent 
or trigger star for mation have be en the subject of debate over the 
past decade (e.g. [McKee & Ostriker^ 2007^ |Vazquez-Semadeni| 
20101 |Hennebelle & ChabrierpOliy The relative importance of 
the various processes is unlikely to be universal: while in some 
regions dynamical effects such as a collision of clouds or conver- 



gence of flows (e.g. |V azquez -Semadeni et al.||2007t [Schneider 
et al.|2QlQ} Duarte-Cabral et al.|2011 ) may be significant, other 



regions can evolve more quietly, from quiescent material possi- 
bly controlled by the local magnetic field (e.g. [Nakamura & Li 
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|2008HAlves et al.||2008 |). The imprints of the initial conditions 
specific to a given cloud are potentially preserved in the prop- 
erties of the gas and dust of their young proto-clusters although 
once star formation begins feedback from the forming stars can 
mask the signatures of these initial conditions. 

Young protostellar outflows may be responsible for clear- 
ing and disrupting protostellar envelopes, limiting the reservoir 
of mass wh ich will end u p being accreted by the central pro- 
tostars (e.g [Fuller & Ladd 2002; A rce & Sargen t 2006; Arce] 
et al.|l2010!). Even though an outflow may carry sufl&cient en- 



ergy to unbind all of an envelope, this energy may well be de- 
posited far from the dense gas, or even outside the molecular 
cloud. Knowing exactly where the outflow momentum and en- 
ergy are deposited is fundamental for understanding whether 
outflows are a source of turbulence capable of providing support 
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against gravity, clearing circumstellar regions and slowing (or 
enhancing) star formation. However, the efficiency with which 
outflows can drive turbulence is strongly debated (se e review 
from McKee & Ostriker 2007 ). The analytical model of Matzner 
I& McKee,(2000j estimates a typical outflow injection scale, the 
scale at which the momentum is most efficiently deposited, of a 
few tenths of a parsec. Numerical models (e.g. Nakamur a & Li 
[20071 [CmoTl et al. 2009, 2010 ) confirm that outflows can drive 
supersonic turbulence at sub-outflow scales and produce smooth 
velocity structures on outflow-scales that sweep up and dissi- 
pate smaller structures, even though outflows are not the primary 
source of turbulence at all scales in molecular clouds. The efl'ect 
of the outflows is to produce an overall flatter density spectrum, 
which limits the fragmentation and infall of material onto the 
final protostars. 

Observationally, the evidence suggests that outflows can 
provide barely enough kinetic energy to support entire clouds 
against gravity (e.g. Arce et al.||2010 |Maury et al.||2009| ), even 
though they may be abl e to maintain the super sonic turbulence 
at sub-parsec scales (e.g. [Nakamura et al.|2011| ). 

Nevertheless, the importance of outflows on the energy 
balance of a region is likely to depend on the physical properties 
of the region. However, it is observationally difficult to measure 
the direct injection of turbulence due to outflows and assess 
its impact on star formation. One reason is that in dynamic 
cluster-forming regions, with large-scale motions, and on-going 
star formation, the direct eflfect of outflows can be easily 
confused and masked by other dynamical events. 

Here we present a study of the B59 star-forming region, 
within the generally quiescent Pipe Nebula. Low-density clouds 
such as the Pipe provide the best opportunities to directly ob- 
serve the impact of outflows where stars have began to form. In 
Sect. [2] we present an overview of the Pipe Nebula and its active 
region, B59. In Sect.[3]we detail the observations, and we present 
the results on the physical properties in Sect.[4j including optical 
depth and estimates of excitation temperatures. The outflows of 
the region are presented in Sect. [5] The ambient cloud's dynam- 
ics and hierarchical structure are described in Sect. [6] (the details 
on the hierarchical structure of B59, performed using a dendro- 
gram technique, is described in detail in Appendix [X|) Finally, 
we discuss the outflow-dense gas interaction in Sect. [7j and we 
outline our conclusions in Sect.[8l 

2. The Pipe Nebula and B59 



Av (mag) 



The Pipe N ebula 
eFaL][j0Q6 



located at ~ 130 pc[^from the Sun j Lombardi 
is part of the Gould Belt of star-forming regions. 
The global features of this cloud appear similar to other Gould 
Belt clouds, with a total mass comparable to the few x 10^ M© 
of the Taurus-Auriga complex, and a filamentary structure with 
a magnetic field threading it (l Alves et"ar] [2008). However the 
Pipe has an extremely low star formation efficiency of less than 
0.1% compared to other Gould Belt clouds which have efficien- 
cies of ~ 2 - 20%. In fact, only one active star f orming clump, 
B59, has been found thus far in t he Pipe cloud ( Brooke et al. 
20071 jForbrich et aL][2QQ9l [2QT0| ). This small group of proto- 
stars within B59 are known to power more than one molecular 



Alves & Franco 



([2007} tance of 145 pc, but 

we adopt 130 pc as it lies within the uncertainties of these authors, and 
it is the distance most commonly adopted in other studies of the Pipe 
Nebula. We note, however, that uncertainties on the distance are of the 
order of 10% and we are therefore possibly underestimating masses by 
20% by assuming the nearer distance of 130 pc. 
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Fig. 1. Extinction map of B59 in grayscale and contours (Roman 



Zuniga et al.| 2009) , showing the known YSOs in the region (from 
Brooke et al.|2007| [Forbrich et al. |2Q09 ), except for the Class III ob- 
jects. Contours range from an Av of 5 to 20 in steps of 5, and then by 
steps of 20. Note that there are only four known Class O and I protostars 
in the region. The black labels show designations of regions within B59 
as adopted throughout the paper. Some cavities are indicated with grey 
arrows and labels, and are discussed further in the text. 



outflow ( ^ishi et a l.'1999 ; Riaz et al.|2009| ). Figure [T] shows the 
extinction map of B59 and identifies several features and regions 
which will be referred to in this work. 

Several scenarios have been suggested for the formation of 
B59 as an active star-forming clump at the edge of the Pipe 
Nebula. One suggests that B59 formed solely governed by grav- 
ity ( Heits ch et al. 2009 ). In another view, the magnetic field plays 
an important role guiding the colla pse of material tow ards the 
main filament ( Alves et al. 2008j). |Peretto et al. (2012) propose 
that the formation of B59 at the edge of the Pipe is due to a com- 
pression front (west of the B59, beyond the western cores) from 
a nearby OB association which has gathered the gas in a "snow- 
plough" fashion. The eff'ects of gravitational focusing playing a 
role later on, with a contraction of the gathered gas and the for- 
mation of stars in the centre of B59. 

Our previous knowledge of the molecular gas structure in 
the Pipe Nebula, however, was limited to the 4' resolution maps 
of the / = 1 ^ transitions of ^^CO, ^^CO and C^^O by 
Onishi et al. ( 1999| ). Beyond this, molecular line observations 
of the Pipe Nebula consist of single spectra of dense gas trac- 
ers such as NH3, CCS, HC5N and other early- and late-time 
molecules (Rathbome et al. 2008t |Frau et al.|[20TQl [2012), and 
C^^O (Muench et al. 2007 ) towa rds the peaks of cores identified 
in r resolution extinction maps ( Lombardi et al.|20Q6||Romai> 
Zuniga et al.p010| ). 



3. Observations 



Using HARP at the JCMT ( [Buckle et al.]2009l ) in May and June 
2010, we have mapped ^^CO and C^^O / = 3 U 2 (at 330.6 and 
329.3 GHz respectively) over the entire B59 star-forming region 
(~ 0.1 1 deg^). The observations were taken with a sky opacity at 
225GHz varying between 0.05 and 0.08. The data were observed 
as several raster maps which have an original spatial resolution 
of 15'' and a spectral resolution of 0.05 km s"^ The r.m.s. noise 
level on the final B59 dataset is 0.22 K (in T* ) with 0.25 km s"^ 
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Fig. 2. Integrated intensity ( J T*j^dv Kkrns"^) maps of the three isotopologues observed. For the C^^O (left) and ^^CO (middle), the integration 
range was 1 to 5.5 kms~^ The ^^CO (right) was integrated from -5 to 15 kms~^ Here we can see that the C^^O follows the shape as seen in 
extinction (Fig.[lJ, while the ^^CO emission starts to pick up the outflowing material, best seen in the ^^CO emission dominated by high velocity 
gas. 
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Fig. 3. Plot of C^^O (purple line), ^^CO (light blue line), with a veloc- 
ity resolution of 0.1 kms~\ and ^^CO (dark blue line), with a velocity 
resolution of 0.3 kms~\ in B59, at the position of the Class source 
Bll, at RA = 17Ml"^23.0% Dec = -27°24'32.8", showing the outflow 
wings and self absorption dip in ^^CO and ^^CO. 



channels. Complementary ^^CO / = 3 ^ 2 (345.8 GHz) data 
were obtained with HARP using JCMT Director's discretionary 
time. These data reached an r.m.s. noise level of 0.2 K (in T^ 
at 0.5 km s"^ velocity resolution. All data were corrected for the 
telescope main beam effi ciency of y/mb = 0-66 ( Buckle et al. 
|2Q09l [Curtis et al.|2010a| ). 

The data reduction was performed using the ORAC-DR 
pipeline ( [Cavanagh et al.|2 008 ), using the recipe reduce science 
NARROWLiNE, with miuor modifications tailored for our dataset, 
with a pixel size of 7.4''. This reduction procedure automati- 
cally fits and corrects baselines and removes bad detectors based 
on the time-series data and on the r.m.s. noise levels of the fi- 
nal maps. The final reduced maps for each area were gridded 
together using Starlink software. 

All datacubes were convolved to a 20'' resolution in order to 
suppress some of the high frequency noise, and were masked to 
pixels where the peak signal to noise (estimated at each pixel) 
was higher than 4 in the smoothed map. This masking procedure 
was particularly eff'ective for the weaker C^^O, but it left virtu- 
ally unchanged the datacubes for the other two isotopologues. 
The reduced integrated intensity maps are presented in Fig. [2j 
and Fig. |3] shows an example of spectra at the position of a pro- 
tostellar source. 



4. Opacity and temperature structure of B59 

4.1. CO-to-H2 relation 

The column density structure of B59 as seen with the dust extinc- 
tion a nd emission is now relatively we ll known ( |Roman-Zuniga 
[eFa l. 2009, 20T2[ [Peretto et al.||2&T2| ). Unlike the dust and gas 
emission, the dust extinction is not sensitive to the temperature. 
As such, studying the relation between the H2 column densities 
derived from the dust extinction and the gas emission can give an 
insight to optical depth eff'ects, temperature eff'ects and the local 
abundance of a molecule with respect to H2. 

Figure |4] presents the relation between H2 column densities 
and the ^^CO and C^^O emission in B59 as scatter plots. The 
H9 column density was deriv ed from the extinction map of B59 
( Roman-Zufiiga et al. 2009) assumin g MH7)Mv = 9.4 x 10^^ 
cm"^ mag-^ (e.g. 'Bohlin et al. 1978] |Pineda et al.||2008] . This 
is plotted against the^CO and C^^O integrated intensities cor- 
rected for the telescope efficiency. Both the extinction and the 
molecular line maps have 20" resolution and were resampled to 
a common area and pixel size, to enable a pixel-by-pixel com- 
parison. 

On the figure, the coloured symbols indicate the correlation 
predicted by the no n-LTE radiative transfer model RADEX ( |van| 
|der Tak et al.||2007) for diff'erent kinetic temperatures, using a 
volume density above the critical densities of ^^CO and C^^O (3 - 
2). The H2 column densities are estimated assuming a fractional 
abundance of ^^CO and C^ ^Q with respect to H 7 of 1.4 x 10"^ 
and 1.7 X 10"^ respectively (^Frerkin g et al.|l982 ). These models 
are only meant to illustrate the behaviour of the emission at dif- 
ferent temperatures, and they do not represent actual fits to the 
datapoints. The curves show the turn-over between a linear cor- 
relation between gas emission and H2 column density and where 
the lines become optically thick. The linear portion of the curve 
is related to the abundance of ^^CO and C^^O relative to H2. 

To estimate the masses and column densities directly from 
the gas emission (Appendix |A]), we estimate a "conversion" 
factor, Xco, on the linear portion of the datapoints as Xco = 
A^(H2)//co (e.g. |Pineda:etaL|2008| , where A^(H2) is the H2 col- 
umn density and Iqo is the integrated intensity of the molecular 
transition. 

For C^^O, a linear fit to all the datapoints provides a Xqi^q = 
1.82X 10^^ cm"^K"^km"^s (purple dashed line, Fig.|4]left panel), 
which will be used to estimate the masses for the bulk of the 
cloud. However, points above and below the purple line, with 
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C180(3-2) Integrated Intensity (K.km/s) 13C0(3-2) Integrated Intensity (K.km/s) 

Fig. 4. Scatter plots of the H2 column density derived from the extinction map against the C^^O (left) and ^^CO (right) integrated intensity 
( J T^b<i V Kkm s~0 for each pixel in the map. The purple dashed-lines represent the linear fit of the data points, with an upper limit constraint on 
the column density for the case of ^^CO. Overlaid on the scatter plots are the predictions of RADEX models at diff'erent gas temperatures (on the 
left: 7.4 K in green and 9.4 K in blue crosses; on the right: 6.6 K in yellow, 9.0 K in green and 11.4 K in blue triangles). For the C^^O (left) the 
green and blue dashed lines represent the linear approximations to the respective colour-coded RADEX model. 
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Fig. 5. Integrated intensity of C^^O (3-2) in the central region of B59 
(gray scale and contours). The protostars are marked with star sym- 
bols, and labeled following [Brooke et al.lpOOT] ). The green and blue 
contoured regions correspond to the regions identified from the C^^O 
scatter plot (see text and Fig.|4]left panel), and the green and blue de- 
partures of the ^^CO integrated emission. Though adjacent, these two 
regions have a temperature diff'erence of ~2 K, the blue region being 
warmer. Note that while the colder green region hosts a single proto- 
star (BIO) and a "starless" C^^O peak, the warmer blue region has three 
embedded protostellar objects (Bll, B7 and B8). 



/cisQ > 1.5 Kkms"\ fall into two spatially distinct regions of 
B59 (Fig. [5] green and blue contoured regions respectively). 
These two regions are of particular interest due to their op- 
tical depth and temperature structure. To better estimate the 
masses and virial parameters in these two specific regions, we 
used a linear fit to the green and blue models, for N(tl2) < 
8 X 10^^ cm"^ and /cisq < 3 Kkms"\ shown as blue and 
green dashed lines on Fig. |4] left panel. These correspond to 



Xci8o(green) = 2.89 x lO^^ cm-^K'^km-^s and Xci8o(blue) = 
1.23 X 10^2 cm-^K-^km-^s. 

For i^CO we estimated Xuco = 2.85 x 10^^ cm'^K'^km-^s 
by constraining the linear fit to pixels where N(tl2) < 3 x 
10^^ cm"^ (purple dashed line on Fig. 4] right panel). We use 
this value to estimate the mass of the bulk of the cloud, but note 
that it will only provide a lower limit in those regions that show 
significant departures from a linear relationship due to optical 
thickening. 



4.2. Optical depth of C^^O and ^^CO 

To measure the optical depth of the ^^CO and C^^O and study its 
variation across the cloud, we have followed the method in Laddl 



et al. ( |1998| ). This uses the ratio of the integrated intensities, and 
assumes the same excitation temperature for both species. The 
method is quite insensitive to the line width assumed, and there- 
fore we have adopted a velocity width of 1 km s"^ throughout. 

The map of the ratio of ^^CO to C^^O integrated intensity 
(not shown) has values ranging from 2 to 4 in the denser regions 
and approaching 9 towards the edges of the cloud. A linear fit 
to the pixel-by-pixel comparison of the two species suggest an 
abundance ratio of ^^CO to C^^O of around 6.4. However this 
value is significantly aflfected by the denser regions where the 
^^CO line becomes optically thick. Therefore, we adopt an abun- 
dance ratio / of ^^CO with respect to C^^O, of / = 8.4 ( Frerking 
et al. 198 2), which is more consistent with the ratio we find at 
the edge of the cloud. This ratio implies C^^O optical depths 
ranging between 0.1 to 0.5 over most of the cloud. However an 
optical depth ranging between 0.5 to 1.25 is found towards the 
central clump (Fig. [S]), consistent with the expected increase of 
opacity from the scatter plot. Therefore, we will consider 
C^^O to be optically thin throughout most of the cloud, but only 
marginally so towards the central cores. 

4.3. Excitation temperature structure 

Since the ^^CO peak emission is optically thick in the bulk of the 
cloud (ti3 > 1), we estimate the excitation temperature of ^^CO 
assuming LTE conditions and optically thick emission. Using the 
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Fig. 6. Left: Colour scale map of the C^^O opacity calculated from the ratio of ^^CO to C^^O emission. The C^^O optical depth is less than 0.5 in 
most of the cloud and only in the central B59 core it rises to higher values, with three local peaks where it reaches values around unity. Centre and 
right: ^^CO excitation temperature map in colour scale. This shows an increase of temperature towards the central B59 (zoomed in on the right 
panel). In all panels, the contours are the dust extinction as in Fig.[T] The darker contour at = 10 which delineates the regions most likely to be 
in LTE. 



radiative transfer equations under this assumption (e.g. |Rybicki| 
|& Lightman|1986| ), we have the relation: 



rex(^'C03-2) = 



15.87 



ln[l + 15.87/(r^,,(i3CO) + 0.045)] 



(1) 



where 15.87 = hvik for the frequency of ^^CO / = 3 ^ 2, 
where h is the Planck constant, v is the frequency and k is the 
Boltzmann constant. Here, T^ax is the peak main beam temper- 
ature of the line. Where ^^CO is self absorbed, the peak used 
corresponds to the absolute maximum of the emission. 

The resulting excitation temperature map is shown in Fig.[6j 
where we can see that most of the cloud is at ~ 9 K, and the star 
forming core is at 10 - 12 K. Given that the ^^CO traces some 
of the outflow emission, in regions with self absorption we will 
be measuring the excitation temperature of the warmer outflow- 
ing gas. As such, outflow shocks can be seen in this image as 
local temperature maxima, even outside the B 59 central region 
(for instance in the SE U-shaped ridge). Outside the contoured 
region shown in Fig. [6] the temperatures decrease to as low as 
5 K. However, in these diflTuse ^^CO emission regions, not only 
is it likely that the ^^CO has lower optical depth, but also the 
densities are likely lower than the critical densities to maintain 
LTE, so the detailed results should be treated with caution. 

However we can compare this temperature structure with 
what is predicted from the ^^CO scatter plots (right panel of 
Fig. |4]). Firstly, we can identify three departures from the lin- 
ear relation, explained by diflTerent gas temperatures (RADEX 
models in yellow, green and blue triangles). The first departure, 
at /i3co = 2 Kkms"\ corresponds to the NE ridge of B59, and 
is consistent with very low temperatures (yellow triangles). The 
RADEX predictions indicate temperatures below 7 K, while the 
LTE estimate (Fig. [6]) indicates temperatures closer to 9 K. 

The other two departures, one where I^co ~ 6 Kkms"^ and 
the other where I^qq has its highest values, at column densities 
below 4x10^^ cm"^, correspond to the same two regions as those 
selected from the C^^O scatter plot (contoured in Fig. 5]). These 
are two adjacent portions in the central B59 region, each contain- 
ing one of the two younger sources in the field (BIO and Bll). 
We find that these two regions show a temperature diflTerence of 
~ 2 K in both molecules (green and blue models on Fig. [4]). Note 



that the warmer (blue) region also includes another two more 
evolved protostars (B7 and B8). This temperature diflference in- 
side the central region is consistent with our LTE ^^CO excita- 
tion temperature estimate shown in Fig.[6j The increase of tem- 
perature in the central region is associated with the position of 
the protostars, and it is likely due to the combined heating eflfect 
from radiative and outflow feedback from these young sources. 
However, the importance of radiative feedback in influencing the 
temperature is only local (within less than 10000 AU from the 
protostars), and even though this could change the fragmentation 
properties of the gas at these scales, it likely does not substan- 
tially alter the fragmentation at the scale of the B59 star forming 
clump (i.e. -0.3 pc). 

5, Outflows 

5.1. Identification of individual outflows 

To identify the population of protostellar outflows in B59, and 
estimate both their global and individual physical properties and 
energetics, we made use of our ^^CO mapping observations. 
These reveal a number of outflows bursting from the central 
sources of B59. The lefthand panels of Fig. [t] show the ^^CO 
blue and red emission from the larger scale outflows in B59. 
Here, we see relatively compact outflow emission bursting to 
the NE of the central region with spatially coincident red and 
blue emission (best seen on the righthand panels of Fig. [7]). This 
emission falls in the extinction Cavity B (Fig.[T]) and the outflow 
appears to trace back to one of the youngest sources in the field, 
the Class O/I source Bll. 

Another set of larger scale outflow features, including a se- 
ries of blue knots, and a red arc can be seen to the SE U-shaped 
ridge. The zoom-in of the central region (right panels) show that 
the driving source of the outflow responsible for the southern 
blue knot seems to be the other young Class O/I source in the 
field, source B 10. On the other hand, the red arc seems to be part 
of the outflow bursting from an older source (source B9). 

This panel also shows that the two Class I sources (B7 and 
B8, light blue stars) sit on a region where the ^^CO from the out- 
flows of the two younger sources (purple stars) start to become 
confused. Therefore, despite the fact that some of the ^^CO wing 
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Fig. 7. ^^CO outflow emission in B59 shown as blue and red contours, together with the YSOs positions plotted on the extinction map from Fig.jl] 
in gray scale (and black contours on the righthand panels). The lefthand panels show the high velocity emission over the entire B59 region, while 
the right panels show a zoom into the central star-forming core. The ^^CO emission was integrated from -5 to 2.7 kms~^ (blue) and from 4.2 to 
15 km s~^ (red). The contours are from 3 K km s~^ on the left panels and from 4 K km s~^ on the right panels, in steps of 2 K km s~^ . 



emission in this particular region is likely to be driven by these 
two Class I sources, it is hard to disentangle from the outflows 
of the other nearby protostars. Finally, the flat spectrum sources 
B1/B2 (the western-most green star) also seem to show some 
signs of outflow emission. Though not very high velocity nor 
very strong emission, these outflow lobes (in particular the red 
lobe) fall near dust Cavity C. 

For many of these outflow features, we cannot clearly iden- 
tify a counter flow lobe. This is likely due to confusion of 
the more compact outflows in the central region. For this rea- 
son, in this paper we will focus on the properties of the un- 
confused outflows as well as on the global outflow properties. 
Interferometry will be important for disentangling these com- 
pact outflows, identifying their driving sources, and provide a 
better estimate of their energetics. 



5.2. Outflow inclination angles 

Before estimating the momentum, momentum flux and energy of 
the outflows, the observed velocities need to be corrected for the 
outflow inclination angle / (where / = 0° is defined to be along 
the line of sight). Observationally, it is hard to infer the inclina- 
tion angle of the outflows with precision, and this introduces a 
critical source of uncertainty, as the velocities (and momentum) 



will be affected by a factor of 1/ cos(/), the energies, by a factor 
1 / cos^(0, and the momentum flux, by a factor sin(/)/ cos^(0. 

Nevertheless, for a few particular outflows in B59, a care- 
ful comparison of the blue-shifted and red-shifted emission can 
shed some light on the overall inclination of the flows. In partic- 
ular. Flow 1 (FigjSj) powered by B 1 1 shows overlapping blue and 
red emission, consistent with being a single outflow lobe, close 
to the plane of the sky. Because the observed maximum veloci- 
ties along the line of sight for the blue and red wings are similar, 
this constrains the outflow axis to being in the plane of the sky 
with less than 5° uncertainty (otherwise the blue and red abso- 
lute velocities would differ by more than a factor 2). In addition, 
the projected opening angle of this outflow is -30°, implying 
that each of the outflow cone walls are inclined by ~15° with 
respect to the plane of the sky, giving / ~ 75°. Note that if the 
outflow axis changes by 5°, the inclination angles would change, 
but while one would increase to 80°, the other would decrease 
to 70°. Because the mass, extent and velocities of the two lobes 
are similar, doing so would increase the estimate of the momen- 
tum and energy by a factor of 3. Therefore, / ~ 75° for Flow 1, 
remains a rather conservative value. 

For Flow 3, the large spatial extent of the collimated blue 
knot (with a faint red counterpart at the end), combined with 
a young driving source (BIO) and a small velocity offset from 
the ambient cloud, also favours an outflow close to the plane 
of the sky. For simplicity, we will adopt the same inclination 
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Table 1. Outflow properties 
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0.63 


0.55 
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Table 2. Momentum flux of individual outflows 







Driving 
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0.19 


-5.5 


75 


0.20 


-21.3 


9.2 


0.15 


0.97 


4.36 


1.1 X 10-^ 


Flow 1 


-Red 


Bll 


0.11 


5.0 


75 


0.11 


19.3 


5.6 


0.14 


0.65 


2.89 


1.2 X 10-4 


Flow 2 


-Red 


B9 


0.23 


5.5 


57* 


0.27 


10.1 


26.1 


0.22 


0.36 


0.46 


1.3 X 10-5 


Flow 3 


-Blue 


BIO 


0.45 


-4.0 


75 


0.47 


15.5 


29.6 


0.17 


0.88 


3.11 


2.9 X 10-5 



* Flow 2 does not show a spatial overlay of the blue and red wings, indicative of it not being as close to the plane of the sky as flows 1 and 3. We therefore adopt a correction factor of 2.9 
for the momentum flux of Flow 2 (i.e. / ~ 57°), based what would be expected from a random distribution of inclination angles jCabrit & Bertout|1992||Bontemps et al.|1996| . If instead 
an inclination of 75° was adopted, the momentum flux of Flow 2 would increase by a factor 5. 
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Fig. 8. Contours of the ^^CO blue-shifted and red-shifted emission (left and right panels respectively). Contours start at 3 Kkms-^ for the blue 
emission and at 4 Kkms-^ for the red, with steps of 2 Kkms-\ overlaid on the extinction map of B59 as from Fig. [l] in gray scale. The black 
arrows show the different individual outflows used to calculate individual outflow properties: Flow 1 (blue and red), Flow 2 (red) and Flow 3 
(blue). 



angle as for Flow 1. Considering the preferential alignment of 
the stronger outflows close to the plane of the sky, we also adopt 
/ ~ 75° to estimate the global outflow properties. 



5.3. Outflow properties 

Estimating the mass of high velocity gas in the outflows, requires 
correcting for the optical depth of the ^^CO linewings, which can 
be done using the ^^CO data. Under LTE and assuming similar 
excitation tempe ratures, a fractional abu ndance of ^^CO relative 
to ^^CO of 62 (L anger&Penzias|l9 93'), and that th e ^^CO wing 
emission is optically thin ( |Cabnt & Bertout||199Q| ), the optical 
depth of the ^^CO (tu) can be determined from the ratio of the 
integrated wing emission of the two isotopologues. 

We integrated the emission below 2.7 kms"^ for the blue 
wing and above 4.2 kms"^ for the red. For each wing, we then 
apply the correction factor of Tu/il - e~'^^^) for the column den- 
sity estimate at each pixel ( Cabrit & Bertout 1990 ; Curtis et al. 
2QlQb| ) . Assuming a kinetic temperature of the outflows of 25 K 
(twice the temperature of the dense gas), a CO abundance with 



respect to H2 of 10 and a distance of 130 pc, the gas mass can 
be written as 



M = 1.78 X lO'^Npi^ < \ Tmhdv > 



(2) 



where M is the gas mass in unit of M© assuming a molecular 
weight of 2.33, Npix is the number of pixels in the outflows, and 
< / Tmbdv > is the average integrated intensity, after the opacity 
correction of the T^b at each pixel. If a higher temperature of 
50 K was adopted for the outflows (as found by, e.g. ,Hatchell| 
|et al.[199 9 ), our mass estimate would increase by 20%. 

To estimate the average momentum, po^i, and kinetic energy, 
Eout, of the outflows, as well as the individual outflow momen- 
tum flux ,Fout. we use 



Pont = m(v)\v - vo\dv 



-f 

^\ J ^(v)(v - 



vofdv. 



Font — Povxhd 



(3) 

(4) 
(5) 
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Fig. 9. C^^O moment maps: mean velocity on the left panel, and velocity dispersion (FWHM) on the right, in colour scale overplotted with the 
integrated intensity in contours, and with the sources plotted as grey stars. The outflows identified in Fig. [8] are also indicated as blue and red 
arrows on the left-hand panel for comparison with the local velocity fields. 



where vq is the velocity of the driving source, m(v) the mass 
corrected for the optical depth, and td is the dynamical time of 
the outflow, defined as = //Vmax, I being the length of the flow, 
and Vmax the maximum velocity of the flow. 

The global momentum and energy of the outflows are sum- 
marised in Table[T] where M is the gas mass contained in each of 
the outflow lobes, /7out is the momentum and £^out is the kinetic 
energy, shown here both without a correction for inclination an- 
gle, and assuming an inclination angle of 75° (pl^ and E^^). 

Table [2] shows the individual momentum flux for the four 
main flows of B59 (Fig [5]). In the table / is the projected length 
of the outflow in the plane of the sky, / is the assumed incli- 
nation angle, Vmax is the line-of-sight maximum velocity with 
respect to the cloud's ambient velocity of 3.5 kms"^ The fol- 
lowing columns give the remaining outflow properties, corrected 
for the inclination angle. The momentum flux of Flow 1 is 
higher than the momentum flux of both Flows 2 and 3 and is 
consistent with the expec ted values for Class sources (i.e. 



10 ^ M© kms yr \ Bontemps et al. 



1996 



Curtis et al. 



M© kms yr ^ Bontemps 



|2010b| ) . Flows 2 and 3 have momentum fluxes simil ar to those 
of a Class I protostar (i.e. ~ 10"" ' ' 
let al.| 19961 [Curtis et al.pOlOb] ). 



6. Ambient cloud 

6.1. Dynamics 

To study the dynamics of the ambient gas of B59, we use the 
combined information from all three CO isotopologues. With 
the C^^O emission in particular, we investigate the velocity gra- 
dients and changes in the velocity dispersion across the re- 
gion (Fig. [9|. The velocity map of the region (left panel of 
Fig. [9]) shows that the NE ridge is uniform in velocity, at around 
3.5 km s"^ Not only do we find little variation in velocity in this 
ridge, but also the linewidth is very small (with a FWHM as low 
as 0.3 kms"\ right panel of Fig.|9]). 

The U-shaped ridge (Fig.[T]) shows a more complex velocity 
structure, with higher velocities than the rest of the cloud (at 
around 4 kms"^). This is somewhat surprising, since this edge 
correlates with a blue outflow knot from Flow 3 (Fig. [7]). This 
strengthens the idea that this is an outflow close to the plane 
of the sky where the redshifted U-shape shows the compressed 



material behind the outflow, although there are no signs of the 
linewidth broadening in this ridge (right panel of Fig. [9]). 

The red arc from Flow 2 is also correlated with a C^^O and 
extinction arc, delineating the righthand side of the U-shape. 
This C^^O gas has higher velocities than the centre of B59, con- 
sistent with being pushed by the red- shifted Flow 2 which has a 
smaller inclination angle than Flows 1 and 3. 

In the central clump of B59 there is a velocity gradient of 
about 1 kms"^ across 0.1 pc (left panel of Fig.|9]). There is also 
a local increase of the linewidth, with the velocity dispersion 
(FWHM) ranging from 0.5 km s"^ to 1.2 km s"^ 

This increase of linewidth can also be seen through position- 
velocity (PV) diagrams (Fig.fTOlandfTT]). The first cut (PV 1) is 



along the axis of Flows 2 and 3. On the left panels of Fig. 
we can see the blue and red outflow knots from Flow 3 at 



11 
l.V 



oflTset as well as the red knot of Flow 2 at -11.5' oflTset, in ^^CO 
and ^^CO. The velocity structure of the U-shaped ridge is seen 
from oflTsets 3' to 14' where, while the eastern end is described 
by velocities ~3.2 kms"^ the U-shaped and western end have 
stronger emission at ~4 km s"^ (cf. Fig. [5]). The ^^CO traces these 
two velocities throughout the ridge, but in fact, the U-shape it- 
self is correlated with the red component (centre-left panel of 
Fig. [To]). Cavity A is therefore also seen as a cavity in the low- 
velocity redshifted ^^CO emission. The central region of B59 is 
delimited by a sharp edge (at an oflTset of 20.5'), where the emis- 
sion from C^^O is broad, and after which the emission drops 
sharply down to the noise in Cavity C. This cavity, seen in ex- 
tinction and with both ^^CO and C^^O, is delimited on both sides 
with high velocity outflowing gas. 



The right column of Fig. 1 1 shows the cut through the NE 
ridge of B59 (PV 2, from 1' to 12' oflTset), confirming the small 
linewidths and constant velocity throughout. The emission from 
the central portion of B 59 is seen at oflTsets between 14' and 20' 
and, similarly to PV 1, is delimited by a sharp edge to the west 
(leading to the same Cavity C as seen in PV 1). The C^^O emis- 
sion shows significant local broadening and the velocities agree 
with those of ^^CO. In particular, the red outflow wings correlate 
with regions where the C^^O is broader with stronger red emis- 
sion. The marked blue outflow wing can be tentatively seen with 
C^^O but it approaches the noise level of the data. Nevertheless, 
the lower velocity blue shifted ^^CO emission at an oflfset of ~20' 
is also seen in the C^^O PV diagram. 
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Dense Gas: CI 80 (3-2) Intermediate densities: 13C0 (3-2) Outflow emission: 1200 (3-2) 
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Fig. 10. Left: Integrated intensity ( J T^^^v) map of C^^O showing the distribution of the dense gas in B59; Centre-left: ^^CO channel map at 
V = 3.7 km s~^ (slightly red-shifted with respect to the cloud's ambient velocities of 3.5 km s~^), tracing the intermediate density gas, and showing 
some of the cavities discussed in the text; Centre-right and Right: ^^CO channel maps at v = 4.8 km s~^ and v = 2.1 km s~^ respectively, illustrating 
the red-shifted and blue-shifted outflows of the region. The position of the two position-velocity (P-V) diagrams of Fig. [TT] are shown as white 
arrows and labeled as PV 1 and PV 2 in the four panels. 



PV 1 TA* (K) PV 2 TA* (K) 




Fig. 11. P-V diagrams of B59 in ^^CO (top panels), ^^CO (centre) and C^^O (lower panels). Each column is the P-V cut as labeled at the top and 
as shown in Fig.[To| The position of some of the outflows are indicated with blue and red dashed lines and cavities are also indicated. The contours 
are from 1 K with steps of 1 K for ^^CO, from 0.3 K with steps of 0.6 K for ^^CO and from 0.25 K with steps of 0.3 K for C^^O (in T^. 



6.2. The hierarchical structure 

To study the hierarchical structure of B 5 9, we have used a den- 
drogram extraction technique developed by Rosolowsky et al. 
( |2008| ). This method is useful to understand the different struc- 
tures that lie within B59. However, like any other clump- 
extraction technique, it cannot be blindly used in a region. In 
the case of B59, most of the leaves retrieved are not "real" cores 
- in the sense of structures which are progenitors of future or 
currently forming protostars. In fact, only two leaves in the cen- 



tral B59 correspond to actual star-forming cores. The remaining 
structures trace the two halves of outflow emission, or less dense 
material which may or may not be bound and form stars in the 
future. A detailed description of the extraction analysis and re- 
sults which confirm the broken structure of B59 can be found in 
Appendix [A| In summary: 

- The SE U-shaped ridge is not a single structure but two dif- 
ferent arch-like structures, consistent with two compressed 
fronts from two diflTerent outflows (Elows 2 and 3). This sug- 
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gests that the gas once in Cavity A has been cleared by Flow 
3 and compressed, forming the eastern side of the U-shaped 
ridge. 

- The Western cores are physically separated from the main 
B59 structure, consistent with the existence of a gas-deprived 
cavity (Cavity C). The presence of outflowing gas possibly 
associated with the B1/B2 system, which is seen at the edges 
of this cavity, suggests that this cavity was possibly cleared 
by the outflow itself earlier in the evolution of these sources. 

- The NE ridge is, at present, a coherent structure, gravitation- 
ally bound, uniform in velocity, and with low velocity dis- 
persion. It could be on the verge of forming future prestellar 
cores; 

- And finally, the dendrogram analysis finds that the central re- 
gion does not show evidence for further fragmentation into 
new cores, apart from the two cores surrounding the younger 
protostars in the field. All the other structures found in this 
central region are components of the outflowing gas, partic- 
ularly along the walls of Flow 1 . 

7. Discussion: Outflow Interactions 

7.1. Outflows and the Cloud 

Our hierarchical study of the region also provides a reliable mass 
estimates throughout the diflTerent scales of the cloud. These can 
be used to estimate the binding energy at diflTerent size scales, 
which can be compared to the kinetic energy carried out by the 
outflows, to estimate their respective impact. The total kinetic 
energy of the outflows (Sect.jsJ is -18 M© km^s"^, assuming an 
inclination angle of 75°. We estimate the binding energy (i.e. the 
potential energy) of the dense material using the C^^O masses. 
For 30 Mq and a radius of 0.30 pc (Appendix [Al), and assuming 
a densi ty profile oc r"^ (similar to that found by ]Roman-Zuniga 
|et al. 2012), the gravitational potential energy is 12 M© km^s"^. 
The outflows therefore carry an energy comparable to the bind- 
ing energy of the dense material. However, the kinetic energy of 
the outflows may be deposited outside the cores, and therefore 
may not be disrupting the dense cores themselves. 

In fact above we have shown that the outflows are correlated 
with the shapes and cavities seen in C^^O and extinction out- 
side the central region, indicating that the outflows are able to 
push, carve and shape the material (e.g. Sect. |6]). For instance 
the material along the SE U-shaped ridge is a fossil structure 
from Flows 2 and 3. Its U- shape, its velocity structure, and its 
local temperature maxima are a result of the impact of outflows 
on an originally less dense medium. 



7.2. Outflows and the Dense Gas 

Despite the clear impact on larger scales of the outflows in the 
B59 cloud, their importance on the dense gas is harder to dis- 
entangle. For instance, in the denser central region of B59, we 
see temperature peaks likely due to the combined eff'ect of ra- 
diation from the individual protostars, and outflow shocks from 
the ensemble of four protostars (Fig. [6]). We detect a broadening 
of the C^^O emission towards the central clump, which could be 
in principle be associated with a combination of rotation, shear, 
infall and outflow motions. But an upper limit for the contribu- 
tion of infall to the linewidth can be estimated by assuming that 
the entire B59 central clump is in free fall collapse. This would 
result on an infall velocity of ~ 0.2 - 0.3 km s"^ (generating a 
velocity dispersion of that order or less), insuflftcient to explain 
the observed line widths. 



To investigate more closely the possible eff'ect of the out- 
flows on the line width of dense gas, we have studied the emis- 
sion in an area where there is initial evidence for the interaction 
of the outflow with the dense gas (in the form of a cavity) and 
where no significant velocity gradients are detected. This is the 
case along Flow 1 , which is an unconfused and well collimated 
outflow located in the densest part of B59 (Fig. IT]). For this test 
case, we compared the average profiles of the integrated 
intensity as a measure of the outflow strength and location; the 
visual extinction as a measure of the quantity of material along 
the line of sight; and the C^^O linewidth as a measure of the 
velocity dispersion. To do this, we took several slices along a 
direction perpendicular to the outflow main axis. 

Since the outflow profile changes along its length, we esti- 
mated the average profiles within two separate regions. These 
are shown as boxes in Fig. [12] (left panels) overplotted on the 
central region of B59 in the integrated intensity (top), the 
extinction map (middle) and the C^^O linewidth (bottom). These 
boxes are colour-coded as the respective profiles, shown in the 
other panels of the same figure. To serve as reference, the dashed 
lines show the position of the peak of ^^CO integrated intensity. 

Closest to the driving source, inside the blue boxed region, 
the outflow is collimated and constrained to a -I'-wide region 



(top row, central panel of Fig. 12). In the second region, de- 
limited by the green box, the outflow is broader (top row, right 
panel). Here, the blue shifted emission starts forming an arched 
lobe opening towards the east, while the the redshifted outflow 
ends with a sharp edge (Fig. [7] and [8|. This sharp tip is immedi- 
ately followed by a linewidth peak, reminiscent of an impacted 
outflow knot (Fig.[T2j first column, lower panel). 

In both boxes, the outflow falls in a similar- sized cavity in 
the extinction (middle panels). This Cavity B (FigJT]) has clearly 
been cleared up by this powerful outflow. Figure 12 (lower row 
panels) shows the impact of the outflow on the dense gas veloc- 
ity dispersion. Even though there is less C^^O gas in this cavity, 
there is still enough material to detect an increase of linewidth 
along the outflow axis. Immediately before and after the out- 
flow axis, towards higher column densities, the linewidths are 
still relatively high (when compared to the outer parts of B59), 
but decrease steadily moving away from the outflowing material. 
This is the case for the regions in the green box between V and 
3', where the C^^O linewidth decreases until it reaches our de- 
tection limit. This suggests that the outflows are having a direct 
impact on both the shape of the region, and the linewidths of the 
denser material in B59. 

An upper limit on the contribution from the outflow- 
generated turbulence in the total velocity dispersion of the re- 
gion, can be estimated by assuming that the initial turbulence 
of the gas in B59 (before any outflow feedback) was as low as 
that currently found outside the central B59 (0.3 km s"^ FWHM 
in the NE ridge, cf. Fig. [9| and Tables [ATT] and [Al2|), and that 
the increase of linewidth along Flow 1 is solely due to outflow- 
injected turbulence. The maximum FWHM detected along Flow 
1 is ~ 1.0kms"\ so we estimate that the outflow-generated tur- 
bulence in B59 is locally of the order of 0.9 km s"^ FWHM (i.e. 
^out _ 0.4 km s"^). This is, however, most certainly dependent 
on the power of the outflow and on the age and mass of the driv- 
ing source. 

Another way to assess the importance of the outflows as tur- 
bulence generators is to compare the tur bulent kinetic ene rgy 
and the outflow energy (similarly to e.g. Arce et aL] 20101) . In 
the central region, the kinetic energy is £'turb = 4 M© km^s"^, 
while the outflow energy, estimated in the same area is E'out = 
0.62 M© km^s"^, before applying any correction factor for out- 
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Fig. 12. Left column: Close up of the central region of B59, where the colour scale and contours are the ^^CO integrated intensity map (top), the 
visual extinction map (centre) and the C^^O line width map (bottom). The regions used to estimate the average profiles (shown in the centre and 
right columns) across Flow 1 are shown as blue and green boxes. The dashed lines show the position of the peak of ^^CO integrated intensity as 
averaged within each box. Centre and Right columns: Average profiles of ^^CO integrated intensity (top row panels), visual extinction (central 
row panels), and C^^O linewidth (FWHM, bottom row panels). Each column corresponds to a diff'erent averaging box, colour-coded as in the left 
column. The dashed lines and gray-shadowed areas mark the peak and extent of the ^^CO profile as seen on the first row panels, respectively. The 
colour-coded shaded areas show the 1-sigma dispersion of each quantity. 



flow inclination. If correcting for a random distribution of an- 
gles and adopting / -57°, this becomes E^^^ = 5.2 M© km^s"^, 
or even E^^^ = 9.2 M© km^s"^ if assuming an inclination an- 
gle of 75° as in Sect. |5] This suggests that the outflows carry 
more than enough energy to provide the turbulent support of the 
central B59. However, since the central region of B59 shows ve- 
locity gradients and is also likely affected by infall motions, it 
is not straightforward to conclude if the observed turbulence is 
predominantly generated by outflows. 

Instead we focus on estimating the fraction of the outflow en- 
ergy deposited in the immediate surroundings of Flow 1 . Taking 
the regions in the outflow walls we estimate from the C^^O that 
the total kinetic energy is ^turb = 0-9 M© km^s"^. From the re- 
spective portion of Flow 1, we estimate that the outflow energy 
is E^^^l = 2 M© km^s"^. As mentioned in Sect. |5]a change in 
the inclination angle adopted for Flow 1 will increase the esti- 
mate of the outflow energy. This together with the fact that the 
total kinetic energy around the outflow does include other terms 



(such as the thermal and isotropic turbulence), suggests that only 
a fraction (less than half) of the energy carried by the outflow 
is deposit ed in the immediate surr oundings consistent with, for 
example, Cunningham et al. ( 2QQ6| ). The remaining energy is de- 
posited outside the star forming core. 

We stress that this is estimate is for the specific case of 
Flow 1 . This eflftciency of outflow energy deposition will likely 
vary depending on the energy of the outflows and the density 
and velocity conditions of the regions with which the outflows 
interact. 

8. Summary and conclusions 

B 59 is a star forming region with a small and young proto-cluster 
at its centre. In its vicinity the gravitationally bound NE ridge ap- 
pears likely to form stars in the future. This ridge has a coherent 
velocity structure, but presently shows very little sub- structure. 
Another region which may eventually form stars are the western 
cores. However the mass reservoir available in this region is not 
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very large and so if it does become active, it will likely form no 
more than a single low mass protostar/protostellar system. 

All the protostars in B59 are located in a C^^O twisted 
filament-like structure (not seen in the dust) in B59's centre. In 
this region, there is a sharp temperature increase of 2K corre- 
lated with the position of the protostars, likely due to radiative 
heating and outflow feedback. Nevertheless, the heating is a very 
localised eflTect, not likely to substantially change the properties 
of the fragmentation of B59 as a whole. 

We identify a number of cavities in the B 59 region that are 
either correlated with outflowing gas, or adjacent to outflow 
knots. Finally, we investigated the relation between the outflow- 
ing gas and the linewidth of the denser gas and found a clear 
correlation with an increase of linewidth of C^^O along the 
outflow axis. 

Our study of B59 suggests that outflows are responsible for 
sweeping up and compressing gas as they make their way out of 
the cloud. We conclude that the outflows are interacting closely 
with the dense material broadening the linewidths of the material 
in the densest gas. Outside of the star-forming core, and as the 
outflows travel through a lower density material, their impact is 
seen as a re- shaping of the gas and dust that surrounds the proto- 
cluster. Their impact on the dense gas is shown to be an eflftcient 
source of turbulence that provides just enough support for the 
cloud against collapse at sub-parsec scales (~ 0.1-0.3 pc). The 
comparable kinetic and potential energies are consistent with the 
idea that the B59 core is long-lived and slowly collapsing. We, 
therefore, c onfirm that outflows are, as speculated by Cove y et aL] 
( 2010 ) and Roman-Zuniga et al."f2012), potentially responsible 
for supporting the B59 clump against collapse. 

A localised line broadening in the dense material is not com- 
monly observed, since C^^O does not trace high velocity out- 
flowing gas. In B59, most of the gas originally had narrow line 
widths, and the turbulence input from outflows is more easily 
detected. In this context, B59 appears in marked contrast with 
regions such as Serpens (e.g. 



Duarte-Cabral et al.l20l"0l ), where 



no clear evidence for the outflows is seen in the C O emission 
Nevertheless, the actual linewidths of C^^O in Serpens (homo- 
geneously around ~ 1-1.5 km s"^ FWHM) are in fact similar to 
those in the central core of B59. Serpens being a much more dy- 
namic region resulting from a collision of clouds (|Duarte-Cabral 
et al. 2011 ) and with a large number of protostars powering 
outflows (more than 20 Class and Class I sources in Serpens 
against four in B 5 9), it may be that we are simply no longer able 
to detect the underlying quiescent gas in Serpens. 

We suggest that outflows are an important source of turbu- 
lence in the gas around star forming cores, in agreement with 
the modeling work from e.g. Matzner^ (2007 j ; Nakamura & Li 
( |2007| ). However, in regions where the initial dynamics of the gas 
implies intrinsically larger linewidths, the contribution from the 
outflow-generated turbulence can potentially become negligible. 
In fact, the final linewidth will only show a measurable broad- 
ening from the outflow contribution (of at least 30%), when the 
injected turbulence is at least as large as the original initial tur- 
bulence of the gas. Assuming that the measured generated turbu- 
lence (cr^^^ -0.4 km s"^) is close to the typical value of Class 
protostellar outflows, the outflow broadening will only be de- 
tectable in regions where the original cr^ of the gas is of the 
order of 0.4 km s"^ (or less). The relative support provided by 
such an input will also vary depending on the mass and the po- 
tential energy available to a given protocluster. In star-forming 
clumps embedded in more massive molecular clouds (i.e. lying 
in a deeper potential well), the gravitational energy will easily 



exceed such outflow-induced turbulent support, and the later will 
not be enough to halt a larger scale collapse. 

We find that the energy carried by the outflows is larger than 
that needed to provide the current turbulent support of B59, and 
that only a fraction, less than half, of the energy of the outflows 
is efl&ciently converted into turbulent motions of the gas in the 
immediate surroundings of the outflows. The remaining energy 
is, therefore, deposited outside the star forming core. In regions 
denser than B59 this efliiciency may increase, as the interaction 
of the outflowing gas and the environment will be more impor- 
tant. On the other hand, the balance of gravitational and turbulent 
energy will also be more important in such denser regions. 
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Appendix A: The hierarchical structure of B59 

A.1. Method: Dendrograms 

To study the cloud structure of B 5 9, we used a dendrogram code 
developed and applied to star forming regions by |Rosolowsky| 
[eFal. (2008 ) that traces the hierarchical structure of the cloud, 
representing the cloud as a tree-like structure. This technique 
begins with the identification of all local maxima in a datacube 
(the leaves of the tree-like structure), and contours down the 
dataset, finding all the underlying structure that connects the 
diff'erent leaves together (i.e. the branches), down to the root. 
Each connecting point between two structures (leaves, branches 
and root) is identified as a node in the tree. The code we used 
to implement this dendrogram technique is fully described in 
[Rosolowsky et al. ( 2008j ). This code also calculates several prop- 
erties of the cloud at each contour level, such as the integrated 
intensity, the size of the structure, the velocity dispersion, the 
mean velocity and mean position. 

We can calculate the masses within the dendrogram code 
by taking a linear relation between the H2 column density and 
the integrated intensity of the molec ular line considered (Xco) , 
following the procedure described in [Rosolowsky et aL]p008| ), 
where 



Table A.l. Properties of the ^^CO dendrogram leaves 



M= 1.84X 10"^^Xco^co 



(A.l) 



where M is the H2 gas mass (in M©) derived from the CO iso- 
topologue in consideration, Xco is the CO-to-H2 conversion fac- 



tor described in § 4.1 (in units of cm ^km ^s), and Leo is the 



luminosity from the considered CO emission (in Kkms"^pc^). 
This luminosity Leo is calculated as L = Fd^, where d is the 
distance to the cloud and F is the integrated flux of the re- 
gion, i.e. the sum of all the emission in the region, calculated 
SisF = 2/ Ti Se^SOy Sv. The value of 1.84 x 10"^^ arises from the 
units conversion, assuming a molecular weight of 2.33. 

The virial parameter a is estimated for each structure at each 
contour level from the virial balance between the kinetic and 
gravitational potential energy, i.e. a = 2Ek/Eg. Assuming a uni- 
form density profile, this translates into 



a = 



MG 



(A.2) 



where CFy is the velocity dispersion and R is the size of the struc- 
ture as estimated by the dendrogram code. Assuming a den- 
sity profile of p oc r"^ would decrease the viral parameter by 
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a factor of 3/5. Without the presence of any additional terms 
in the virial equation (such as magnetic field or external pres- 
sure), virial equilibrium corresponds to a = 1. However, given 
the uncertainties associated with these assumptions, we consider 
a structure to be gravitationally bound when the virial parameter 
is less than 2. The use of such a parameter is best seen as an in- 
ternal means of comparison of structures within the cloud rather 
than an absolute measure of the individual equilibrium states. 

We performed this dendrogram analysis of B59 on the 3D 
datacubes of the intermediate density tracers ^^CO and C^^O. 
We considered the minimum height of a leaf to be one-sigma of 
the r.m.s. noise level, and chose t he bijection method to ca lculate 
the cloud properties (similarly to [Rosolowsky et al.|2008 ). 



A.2. Results 

A.2.1 . The ^^00 and C^^O dendrograms 



The dendrograms for ^^CO and C^^O are shown in Fig. A.l left 
and right panels respectively. This hierarchical stratification of 
^^CO finds the entire structure of B59 to be interconnected and 
part of the same cloud, while with C^^O the cloud is no longer 
interconnected (likely due to signal to noise limitations). Despite 
the non-interconnectivity, the C^^O hierarchical stratification is 
consistent with that found with ^^CO though with less branching. 
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Fig. A.l. Dendrogram of the ^^CO (left) and C^^O emission (right) in B59, with the virial parameter calculated for each contour. The regions 
where no virial parameter is shown (e.g., the tips of the leaves) are due to the lack of sufficient quality data to calculate the properties, such as 
insufficient number of pixels. The numbers refer to nodes and leaves as found by the code. 



Th e pro perti es of the dendrogram leaves are summarised in 
Tables |A.1| and |A.2| The first column indicates the region to 



which each leaf belongs and the second column shows the den- 
drogram ID number (cf. Fig. A.l A.2 and |A.3|). The follow- 



ing columns show the properties derived at the lowest contour 
of each detected leaf: R is the radius; vq is the mean velocity, 
weighted by the intensity (with an uncertainty of 0.1 km s"^); CTy 
is the velocity dispersion; M is the leaf mass (eq. A.l ); and a 
is the virial parameter (eq. |A.2| ). Bear in mind that the virial pa- 
rameters shown in this table are estimated at the bottom of each 
leaf, but they evolve throughout the hei ght of the leaves, gen- 
erally approaching sub virial values (Fig. |A.1| ). In addition, the 
row before-last shows the properties of the branch that covers 
the entire central region, and the last row shows the total mass, 
equivalent size and virial parameter of the entire B59. For ^^CO, 
the properties of the entire B 59 were measured at the root of the 



dendrogram (root 44, Fig. 



A.l 



left panel), while for C^^O this 
le information of the three dis- 



was estimated by combining t 
connected roots (Fig. |A.l [ right panel, roots 50, 47 and leaf 3). 

These estimates do not account for the optical depth, which 
will particularly aff'ect the central region of B59. Figure|4]quanti- 
fies the possible eff'ect of optical depth for C^^O as deviations at 
higher column densities above the green dashed line. This would 
account for a diff'erence of ~ 15% on average in the column den- 
sity in this central region, which translates into an equivalent 
increase of the mass estimate (and decrease of a). 

A.2. 2. The leaves of B59 

The densest gas within the central clump of B59, is solely de- 
scribed by two leafs in both molecular lines, leaves 4 and 6, both 
gravitationally bound, and both at the systemic velocity of the 
cloud (in C^^O). These correspond to the two cores surrounding 
the youngest protostars of the region: B7, B8 and B 1 1 within leaf 
4, and BIO within leaf 6. Apart from these, we have four leaves 
with ^^CO in the c entral clump (leaves 7, 8, 17 and 18, labeled 
in red in Fig. A.2 ). These are a result of the ^^CO self absorp- 
tion and outflow wings, which produce an overly sub- structured 
tree tracing outflowing gas at opposing velocities, and artificially 
separated by the self-absorption dip. These leaves do not have 



the systemic velocity of 3.5 kms ^ but they vary from 2.9 to 
3.9 km s"^ and trace the dense outflowing gas from Flow 1. 

The SE U-shaped ridge is divided into two halves in both iso- 
topologues (labeled in yellow in Fig. |A.2| and |A.3| ). The lower- 
left half of the U-shaped ridge is identified with leaf 3, at slightly 
redshifted velocities, consistent with the velocity map of C^^O 
(Fig. [9]). The other half of the ridge is found as a red-shifted leaf 
in ^^CO (leaf 5), with a high virial parameter, consistent with 
compressed gas from the redshifted outflow (Flow 2). 

The entire NE ridge is found as leaf 13 (labeled in green in 
Fig. |A.2| and |A.3| ), without any sub- structure within it, and at the 
systemic velocity of the cloud (-3.5 km s"^). This is the only leaf 
in B59 that appears with a a below 1 in both ^^CO and C^^O. 

On the other hand, the western cores branch out from the 
main tree and form a sub-structure of their own. They include 
leaves 10 and 15 (labeled in blue in Fig. A.2 and A.3 ), detected 
with both molecules. The mass and virial parameter of leaf 15, 
in particular, indicates that it could collapse to form future pro- 
tostars, but at present it presents no evidence for ongoing star 
formation. Despite having a similar mass to the central cores of 
B59, its larger size (and therefore lower densities) are indica- 
tive of a less compact core, which has at its disposal a smaller 
mass reservoir from where it could accrete. If it does become 
active, it will likely form no more than a single low mass pro- 
tostar/protostellar system. Spatially nearby, there is the blue- 
shifted ^^CO leaf 16, and the C^^O leaf 21 , both with a high virial 
parameter. These two leaves lay at the edges of Cavity C, sug- 
gesting that they are part of the material pushed out of the cavity 
by an outflow. This gas-deprived cavity is now what physically 
separates the central region and the western cores. 

Finally, we hypot hesise that the ^^CO leaves 9, 1 1 and 14 (la- 
beled in purple in Fig. |A.2[ ) are counterparts of the main outflows 
of the region, but since they are not at very large velocity off'sets 
and since we cannot trace the material between these knots and 
the respective driving source, it is not straightforward to under- 
stand their origin. Leaf 19, in particular, appears isolated and is 
a possible outflow bullet, with the highest virial parameter and 
the highest blue shifted velocity in the cloud (at -2.2 km s"^). 
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Fig. A.2. Six channel maps of the ^^CO emission, T^, whose velocities are shown at the top-left corner of each panel. The plotted grayscale is 
that used for all panels with the exception of the first (top-left panel), where it is from 0.1 to 1.5 Kkm s~^ The black contours show the positions 
of the leaves from Fig. |A.l| whose numbering is colour-coded to visually separate the distinct areas: central region in red, NE ridge in green, SE 
U-shaped ridge in yellow, western cores in blue, and loose knots in purple. 
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Fig. A.3. Four channel maps of the C^^O emission in T^ (in gray scale), whose velocities are shown at the top-left corner of each panel. The black 
contours show the positions of the leaves from Fig.|A.l| and are numbered accordingly. The colour code is the same as for Fig.|A.2| 



A. 2. 3. The masses 

A brief comparison of the masses derived from the ^^CO and 
C^^O emission shows that the structures which are common have 
similar mass estimates. C^^O recovers more mass from the NE 
ridge and from some of the cores, as expected from the higher 
optical depth of ^^CO. Nevertheless, both lines show a total of 
~ 10 Mq in the leaves. However, they differ in the total mass of 



the cluster by a factor of -2.5 (~ 74 M© in ^^CO versus ~ 32 M© 
in C^^O). Such a difference is mostly due to the diffuse gas traced 
only by ^^CO. That said, the fraction of mass in the leaves is 
-11 % using the ^^CO, and ~ 33% for the C^^O estimate. 

Considering only the central clump of B59, our best mass 
estimate is retrieved from C^^O, as estimated at the lower level 
of branch 45. This central region is found to have a mass of 
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-16 M© within a region of ~ 0.1 pc radius, with an average 
velocity dispersion (ctv) of 0.3 kms"^ (i.e. a = 0.7). 

When using the momentum maps constructed over the entire 
velocity range, we measure within a central region of ~0.1 pc 
radii a mass of 10 M©, a mean velocity of 3.5 km s"^ and a mean 
(Ty of 0.35 kms"^ Assuming the uniform profile as that taken 
for the dendrogram analysis, we calculate a virial parameter of 
1.5. The smaller mass estimated using the total integrated in- 
tensity is not surprising since we are now taking a circular area 
around the central region of 0.1 pc radii, out of which not all 
pixels are included in the calculation due to signal to noise lim- 
itations. For the same reason, the higher velocity dispersion is 
also compatible with the dendrograms, likely arising from the 
fact that the larger line widths in B59 are concentrated toward 
the central region. This is also clear lo oking at the variations of 
the virial parameter of branch 45 (Fig. A.2 right), which starts 
at ~2 at higher levels, decreasing down to 0.7 at the lower level 
of the branch. 
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